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I ABSTRACT 
i 

1 We use a representative sample of 65 galaxy clusters observed with the Einstein 

2 IPC to constrain the range of cluster X-ray morphologies. We develop and apply 
" quantitative and reproducible measures to constrain the intrinsic distributions of (i) 

emission weighted centroid variation w^, (ii) emission weighted axial ratio ?y, (iii) 
emission weighted orientation ^o, and (iv) measures of the radial fall-off, a and fi. 
Eor each cluster we use a Monte-Carlo procedure to determine the effects of Poisson 
noise, detector imperfections, and foreground/background X-ray point sources. 

We then use the range of cluster X-ray morphologies to constrain three generic 
cosmological models (0=1, Oo=0.2, and Oo=0.2 & Ao=0.8). Eor each of these mod- 
els, we evolve eight sets of Gaussian random initial conditions consistent with an 
effective power spectrum P{k) oc on cluster scales. Using this sample of 24 
numerical cluster simulations which include gravity and gas physics (but no cooling 
or ejection from galaxies), we compare the X-ray morphologies of the observed and 
simulated clusters. Specifically, we build artificial ensembles with the same distribu- 
tions in the number of cluster photons. X-ray temperature, and cluster redshift as 
the Einstein ensemble; we then compare the observed and simulated distributions in 
, ?7, and a. 

The comparisons indicate that: (i) these three morphological characteristics are 
sensitive to the underlying cosmological model, and (ii) galaxy clusters with the 
observed range of X-ray morphologies are very unlikely in low Oq cosmologies. The 
analysis favors the 0=1 model, though some discrepancies remain. We discuss the 
effects of changing the initial conditions and of including additional physics in the 
simulations. 

Subject headings: cosmology: theory — galaxies: clustering — hydrodynamics — 
intergalactic medium — methods: numerical — X-rays: galaxies 



I. INTRODUCTION 

Over the past decade, studies have provided ample evidence that a significant fraction 
of galaxy clusters have undergone recent growth (Geller & Beers 1982, Dressier & Sheet- 
man 1988, Eorman & Jones 1990). Analytic work by Richstone, Loeb, & Turner (1992) 
demonstrates that the present epoch dynamical state of galaxy clusters is sensitive to the 
cosmological density parameter Oq. We have previously shown that the structure of the 
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X-ray emitting gas in clusters reflects this dependence (Evrard et al. 1993). Thus cluster 
X-ray morphologies provide an important cosmological constraint. 

Here we use quantitative and reproducible measures of cluster X-ray morphology for a 
representative sample of galaxy clusters to evaluate these cosmological constraints. We flrst 
develop and test measures of the cluster centroid variation (Mohr, Fabricant, & Geller 
1993; hereafter MFG93), the axial ratio ?y, the orientation ^o, and the measures of radial 
fall-off, a and fi. We apply these measures to a representative sample of cluster X-ray 
observations to determine the intrinsic cluster distribution in these quantities (§2). Finally, 
we directly compare the morphologies of observed clusters to a set of clusters simulated within 
three different cosmological models; we use the results to place cosmological constraints 
(§3). §4 contains a summary and further discussion of the results. (A Hubble parameter of 
Hq = 50 h'^f^ km s~"^ Mpc~"^ is assumed throughout.) 

2. THF MORPHOLOGIFS OF OBSFRVFD CLUSTFRS 

This section is a description of the measures of cluster X-ray morphology. We include a 
discussion of the observational sample taken from the Einstein archive, and a brief review of 
the standard IPC reduction procedure. Case by case, we provide the measurement details of 
four morphological characteristics: (i) the centroid variation w^^ (ii) the axial ratio ?y, (iii) 
the orientation Og^ and (iv) two measures of the radial fall-off, a and /3. We review the tests 
of our measurement accuracy and display the results. 




Figure 1: Completeness of the Emsiem sample. Col- 
umn 1 contains distributions of 2-10 keV X-ray lu- 
minosity for (bottom) the 55 members of the flux 
limited Edge ei al. (1990) sample and (top) the 
58 Emsiem clusters with published luminosities (Edge 
ei al. 1990, David ei al. 1992). The dashed lines 
in the upper plot form the distribution for the 41 of 
58 Emsiem clusters contained in the Edge ei al. flux 
limited sample. Column 2 contains an analogous ar- 
rangement of the distributions in X-ray temperature. 
The KS test demonstrates that there is a reasonable 
probability that these three samples are drawn from 
the same parent distribution. 



LoglO(Lx{3_io))[erg s ^] Tx[lceV] 



2.1 THF CLUSTFR SAMPLF 

Cluster X-ray morphologies are quite varied. To probe the full range of morphologies it 
is imperative to have a large, representative sample of high signal-to-noise cluster images. 
In the near future the ROSAT public archive will be a valuable resource in this regard. 
For the present study we extract and reduce Einstein Imaging Proportional Counter (IPC) 
images of 65 galaxy clusters. We choose these particular clusters in an attempt to construct 
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a representative sample and because they are among the best cluster images obtained with 
the IPC. Our ensemble of observations contains images with from 10^ to 2 X 10^ cluster 
photons; the median number of cluster photons per image is ~5000. Because this cluster 
sample is not manifestly complete with respect to volume or flux limits, we evaluate the 
effects of incompleteness by comparing the sample to the X-ray flux limited sample of 55 
clusters constructed by Edge et al. (1990). 

Of our 65 clusters, 58 have published 2-10 keV X-ray luminosities {Lx) and X-ray tem- 
peratures {Tx) (Edge et al. 1990; David et al. 1992). 41 of these 58 clusters are in the Edge et 
al. sample. Eigure I compares the distributions in Lx (left column) and Tx (right column) 
of our ensemble (top) with the Edge sample (bottom). A first glance at the distributions 
may lead to the conclusion that the entire Einstein sample has more low luminosity and low 
temperature clusters than both the Edge sample and the Einstein subsample. However, the 
KS test demonstrates that the differences in these distributions are all consistent with sta- 
tistical fiuctuations. Eor each combination we include the KS parameter D and the implied 
probability that the two distributions are consistent. Eor Tx- Edge/Ein58 D=0.I6 (45%), 
Edge/Ein4I D=O.IO (98%), Ein4I/Ein58 D=0.I5 (69%). Eor Lx: Edge/Ein58 D=0.2I 
(15%), Edge/Ein4I D=0.I4 (71%), Ein4I/Ein58 D=0.I5 (67%). On the basis of these re- 
sults, the Lx and Tx distributions for both the entire Einstein sample of 58 clusters and 
the Einstein subsample of 41 clusters are indistinguishable from an X-ray fiux limited sam- 
ple. 

Naturally, we would like to have images of comparable quality and depth for all 65 
clusters. However, variations in exposure time and surface brightness within our sample 
complicate our goal of extracting the intrinsic range of cluster morphologies. Using mor- 
phological measures which are most sensitive to the bright cluster cores reduces sensitivity 
to these sample variations. Each of our measures (except for the radial fall-off) is emission 
weighted and therefore most sensitive to the highest signal-to-noise regions of the X-ray 
images. Our measures are relatively insensitive to the inclusion or exclusion of outer cluster 
regions. Indeed, for a cluster with the mean radial behavior (/3 = 0.65), the emission weight- 
ing scheme outside the cluster core is essentially a radial weighting by where R is the 
projected distance from the peak in the X-ray surface brightness. 

2.2 IMAGE REDUCTION 

We filter the observation to produce a 0.3 to 3.5 keV X-ray image blocked from the 
original 8"pixels to I6"pixels. We then subtract an approximation to the image background 
using a scaled Deep Survey background map; the background is scaled using either the 
observation livetime or the X-ray fiux in regions of the image outside the cluster. Then we 
apply corrections for telescope vignetting and small scale IPC non-linearities (see Soltan & 
Eabricant (1990) and MEG93 for a discussion of the small scale correction). Einally, we block 
the image to 32"pixels and smooth with a Gaussian of 1.5' EWHM. The resulting image has 
a spatial resolution of 2.4' EWHM. Our analysis of these images is confined to the detector 
region within the IPC ribs. 
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2.3 CENTROID VARIATION 

We analyze each X-ray image in annuh of constant width (1.07') and increasing radius; 
as described below, these annuli are not concentric. Because an a priori choice of image 
center biases the results and the cluster centroid often varies as a function of scale, the 
first step in the analysis is defining a cluster center for each annulus. We use a simplex 
minimization routine (Press et al. 1988) to find the annulus position which minimizes the 
distance between the annulus centroid x and the geometric center of the annulus. x is the 
intensity weighted first moment of the annulus 



^^n^x^, -J^^n^y^^ where N = 



(2.1] 



Hi is the number of photons in pixel i and the summation occurs over all pixels within the 
annulus. This minimization is roughly analogous to balancing a donut on an anthill- the 
final position of the donut is the position which brings the plane defined by the donut as 
nearly parallel to the ground as the anthill allows (see MFG93 for more details). Once the 
appropriate annulus position is determined, we use the centroid x calculated at that position 
as the image centroid for that annulus radius. 

The annulus radius is increased and the identical prescription followed until measure- 
ments cover the useful portion of the X-ray image (determined by the IPC ribs and any 
bright point sources). Thus we extract a radial function in the centroid x. To examine the 
variation in this function while minimizing the noise introduced by the low signal-to-noise 
region of the image, we calculate an emission weighted centroid variation over the region 
of the image which satisfies a mean signal-to-noise constraint ((■^) > 5 where <(-^)> is the 
mean signal to noise of all the pixels within the annulus). Specifically, the emission weighted 
centroid variation is 



<=\1^^A 2^N,{x,-{x)y (2.2) 
J / J 

where the summation occurs over all annuli with <(-^)> > 5, Nj is the number of photons 
within the annulus j, and (x) is the emission weighted average centroid of the image. 

We take a Monte-Carlo approach to determine the significance of the centroid variation. 
We ask, "Is the centroid variation we detect consistent with variations introduced by Poisson 
noise, foreground/background point sources and detector imperfections?" To answer this 
question, we first fit our cluster image to a symmetric /3 model of the form 

(I) ,,3, 

where /q, i?c, V P ^^'^ the normalization, core radius, axial ratio and radial fall-off. 
Two other free parameters, a cluster orientation angle and emission center on the IPC, are 
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included because the detector response varies over its face. This /3 model has no intrinsic 
centroid variation. 

We then simulate an observation of this model, including X-ray background, Poisson 
noise, detector response variations and randomly positioned point sources consistent with 
the EMSS log A^-log S distribution (Gioia et al. 1990; see MFG93 for details). We reduce the 
simulated image ("/3-image" hereafter) just as we reduce the original image and carry out the 
measurements described above. By reimaging the fi model from 200 to 1000 times we build 
up a sample of /3-images with a range of morphologies (and centroid variations) due solely 
to non-cluster effects. We use this distribution of centroid variations in two different ways: 

(i) we use the mean of the distribution {w^} j^^q to correct the centroid variation measured 
in the original image for the Poisson, point source and instrumental contributions {w'i = 
^meas ~ {^"^) MC ^^^^^^ ^meas ^1 ^he Original and corrected centroid variations) and 

(ii) we use the width of the Monte-Carlo distribution as one component of the uncertainty 
in the centroid variation. 

The value of depends on the region of the cluster used to calculate it (as discussed 
later, the orientation, axial ratio, and radial fall-off share this characteristic). This depen- 
dence can be produced by a cluster cooling flow located off center with respect to the large 
scale structure of the cluster. In this case the cluster centroid would vary significantly near 
the cluster core and then vary little outside the core. Thus, would decrease as we included 
more and more annuli from the outer regions of the cluster. Because the centroid variation 
depends on the radial scale of the measurement, the choice of an appropriate scale is cru- 
cial. Here we interpret a centroid variation as an indicator of dynamical youth (MFG93, 
Mohr & Evrard, 1994). Thus, we take the radius where reaches maximum significance 
[max{w^/crw}^ where (Ty, is the centroid variation uncertainty). 

In §3 we demonstrate that clusters simulated with an Einstein-deSitter cosmological 
model have morphologies very similar to observed clusters. Thus these cluster simulations 
provide an ideal medium to test our measurement accuracy. We compare measurements 
made on perfect images of the simulated clusters (no instrumental effects or Poisson noise) 
with the measurements made on artificial IPC images of these clusters (see §3 for details 
about these artificial images). In particular, we use a sample of more than 200 cluster images 
with distributions in signal-to-noise and spatial scale similar to those of the Einstein IPC 
sample. 

Eigure 2 contains the results of these comparisons. The first column contains a measure of 
the accuracy of w^; the top row contains a plot of the fractional error as a function of the real 
value. Eor a centroid variation w^< 0.3', the accuracy decreases because of the resolution and 
noise characteristics of our image ensemble (and because we plot the fractional error which 
is proportional to w^~^). The large number of photons in a typical image (~5000) allows 
us to extract variations which are ~ |- the resolution scale of the images. Greater resolution 
would require either IPC images with significantly more cluster photons or an instrument 
with improved imaging and background characteristics (e.g. the ROSAT PSPC). The bias for 
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Centroid Variation Axial Ratio Orientation Figure 2: Measurement 

accuracy evaluated using 
a set of over 200 artificial 
IPC images. The first 
row contains (from left to 
right) (i) the fractional 
error in the centroid 
variation (Aw/wreai = 

(ws - WreaO/Wreal) Ver- 

sus the real centroid vari- 
ation Wreai, (ii) the frac- 
tional error in the ax- 
ial ratio (Ari/rireai = 

(^^meas eal) / eal) Ver- 

sus the measured axial 
ratio rjmeas, and (iii) the 
absolute error in the ori- 
entation angle (A6 in de- 
grees) versus the mea- 
sured axial ratio rjmeas- 
The mean errors for 

-5 5 -2 2 -5 5 

the three parameters are 

(Aw/wreai) = -7.6% (for all Wreal > 0.3'), (Arj/rjreai) = "0.7%, and {A0) = -0.85°. The second row 
contains the histogram of scaled errors (i> = A/ a, where A is the error and a is the uncertainty) for each 
measured parameter. The best Gaussian fits to the i> distributions have the widths a^, = 1.34, cr^ = 0.73, 
and ag = 1.49. For each parameter we use these widths to correct our uncertainties in order to obtain a i> 
distribution of nominal width. 




images with Wj.^al > 0.3' is {Aw / Wj.^al) = —7.6%; on average we tend to underestimate the 
magnitude of the centroid variation. This bias results from our method of correcting for the 
Poisson and instrumental contributions to Wmeas- As described above, we postulate that the 
measured centroid variation has an intrinsic and a noise component: tw^gas ~ "I" (^f )mC" 
This approximation works well at large centroid variation but leads to an underestimate 
for small variations. Clearly, Poisson noise and instrumental effects can also mask a cluster 
centroid variation! 

The bottom row of the first column of Figure 2 contains a scaled error distribution for 
all images [p = A/a where A is the measurement error and a is the measurement uncer- 
tainty). We obtain this distribution by including three contributions in the centroid variation 
uncertainty: (i) the width of the Monte-Carlo distribution of w^^ (ii) the formal uncertain- 
ties which follow directly from Equation 2.1 and Equation 2.2, and (iii) a fO% fractional 
uncertainty which accounts for the uncertainties in the IPC response. This fiat 10% frac- 
tional uncertainty removes the outliers in our scaled error distribution. A Gaussian fit to 
the p distribution demonstrates that we underestimate the uncertainties by 34% (the best 
fit cr = 1.34). We increase all uncertainties by 34% to correct for this underestimate. 

Unfortunately, it is not possible to extract a centroid variation from all 65 cluster images. 
Table I lists for 46 clusters and the maximum centroid shift Ax within the ("^) > 5 region 
for 58 clusters. The table contains the centroid variation and uncertainty in arcmin (tWf ['] 
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Centroid Variation 





Axial Ratio Radial Fall off Figure 3: X ray 

morphologies of observed 
clusters. From left to 
right is (i) the distribu- 
tion of emission weighted 
centroid variation in 
kpc for a sample of 43 
clusters, (ii) the distribu- 
tion of emission weighted 
axial ratios rj for a sam- 

50 100 150 200 250 4 -6 -8 1 A 6 8 1 12 1-4 pie of 5 1 clusters aud (111) 

w^[kp(-] 77 ,6 ^Yie distribution of radial 

fall-off /3 for 48 clusters. The mean and root mean square width of the three distributions are: a) {wg) = 41.3 
kpc and RMS„ = 47.4 kpc, b) (ry) = 0.80 and RMS„ = 0.12, and c) = 0.65 and RMS/j = 0.16. The 
dashed line histogram in the rj distribution contains 38 clusters. The 13 clusters removed from the original 
sample have central cooling times significantly less than 10"'^'^ years (Edge, Stewart, & Fabian 1992, Stewart 
ei al. 1984). 

and (Tw)^ the centroid variation and uncertainty in /i^Q^kpc {w^[kpc\ and cr^^), the radius (in 
arcmin) of the maximally significant wy^ measurement (i?[']), and the maximum centroid shift 
(Ax). Clusters with a "-" entry have centroid variations consistent with zero (where {w^} j^^q 
is greater than ti'^eas)- ^ "*" marks those cluster where the {jj-} > ^ region is too small 
to make a meaningful measurement, and a "p" marks clusters with bright point sources 
near their emission peaks, making difficult to measure with the Einstein IPC. (The higher 
resolution ROSAT PSPC images of these clusters may allow accurate measurements after 
point source removal.) 43 of these 46 clusters have measured Lx and Tx- A KS test 
demonstrates that this sample of 43 clusters is statistically indistinguishable from the fiux 
limited sample of Edge et al. [Tx: D=0.20 (28%), Lx: D=0.23 (15%)). Figure 3 shows the 
distribution of w^. The mean and root mean square width of this distribution are (w^) = 41.3 
kpc and RMS«, = 47.4 kpc. Notes on individual measurements are in Appendix A. 

Because the centroid variation is a signature of dynamical youth, we can use these 
measurements to examine the fraction of substructure in galaxy cluster cores. Of the 46 
clusters with measured w^^ 71% have centroid variations which are at least 3(7^, from zero; 
48% have centroids variations which are at least bcr^, from zero. The Monte-Carlo method 
we use to determine the uncertainties in our measurements also gives us the probability 
that the measured centroid variation in a cluster is merely due to Poisson noise, detector 
imperfections, and X-ray point sources. The probability that is due to an intrinsic, 
cluster centroid variation is > 99.9% in 46% of the sample and is > 99% in 61% of the 
sample. Thus, we conclude that the IPC images of these 46 clusters provide evidence for 
core substructure in from 50% to 70% of the sample, depending on the confidence criterion 
adopted. 

We include Ax because this scale is readily apparent in a contour plot of the cluster 
X-ray emission. Table 1 lists the values for all clusters with at least two meaningful centroid 
measurements within the {jf} > 5 region of the image. We avoid using Ax for comparison of 
X-ray morphologies because it is a noisier measure than (dependent on only two centroid 
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Table 1. Centroid Variation 
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measurements rather than an ensemble of >4 measurements). 

2.4 AXIAL RATIO 

The methods of measuring the cluster axial ratio and the centroid variation are similar. 
Given the correct position for an annulus, we expand the photon distribution defined by the 
annulus in a finite Fourier series. 

4 

/(6)) = ^ A„ cos n6' + 5„ sin n6' (2.4) 

ra=0 

The series is truncated at n = 4 to reach a balance between execution speed and accuracy 
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in the low n terms. The axial ratio comes from the coefficients of the 29 terms: 



^_ ^ r minor ^ r_ oi_ ^ ^ _ V ^ ^ ^^ 2.5) 

rmajor r + Sr 



A formal uncertainty for the axial ratio follows from its dependence on A2 and B2 and 
the appropriate diagonal elements of the covariance matrix from the expansion. The final 
uncertainty for the axial ratio at a particular radius is the sum in quadrature of the formal 
uncertainty and a more realistic uncertainty based on the Monte-Carlo approach (§2.3). 
The Monte-Carlo uncertainty is the RMS variation of the axial ratio over the ensemble of 
/3-images. 

With this approach we measure the cluster axial ratio (and uncertainty) as a function of 
radius. We use this function within the ("^) > 5 region of the cluster image to calculate r], 
the emission weighted average axial ratio. 



E 



N 



.2 



(2.6) 



where N j is the number of photons scaled by the axial ratio uncertainty in each annulus 
and rjf is the measured axial ratio in each annulus. As for t/;^, we use a large ensemble of 
numerical cluster simulations to evaluate the accuracy of our measurements. The second 
column of Figure 2 contains the results for rj. The top row contains a plot of the fractional 
error in ?y as a function of the measured axial ratio. There is no significant bias in this 
measurement; the mean error in the more than 200 images is {Arj /r]j.^al) = —0.2%. The 
width of the best fit Gaussian to the scaled error distribution in the bottom row is (7=0.73, 
27% less than the nominal width. We reduce the uncertainties by 27% to correct for this 
overestimate. 

Columns labeled rj and cr^ of Table 2 contain the list of axial ratios and uncertainties 
for 57 clusters. The columns labeled R['] contain the radius (in arcmin) within which we 
calculate the emission weighted average. 51 of these clusters have measured Lx and Tx 
which we use to compare this sample to the fiux limited sample of Edge et al. The KS test 
demonstrates that this sample of 51 clusters is statistically indistinguishable from the Edge et 
al. sample {Tx'- D=0.19 (31%), Lx'- D=0.24 (10%)). Eigure 3 contains a histogram of rj for 
those 51 clusters; the mean and width of this distribution are = 0.80 and RMS,, = 0.12. 
Appendix A contains a detailed discussion of individual measurements. 

McMillan, Kowalski, & Ulmer (1989- hereafter MKU) study the orientation and axial 
ratio of the X-ray emission from 49 Abell clusters. They focus on the faint outer region of 
the cluster image. Specifically, the region they use is defined by two isophotes: a low fiux 
level just above the noise and a high fiux level defined so the entire region contains 20% of 
the cluster fiux. Because we emphasize the core of the X-ray emission, our results differ 
from those of MKU in cases where the morphology of the X-ray emission varies with radius. 
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Table 2. Radial Fall-off, Axial Ratio and Orientation 
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central cooling time < 10^'^years * member of Edge et al. sample 

+ /3's accurate to 17% "* no measured X-ray temperature 



The measurements on three clusters make the differences abundantly clear: (i) the X-ray 
emission in Abell 754 is very flattened near the core and becomes more circular outside the 
core [rj = 0.579 versus rjMKU = 0.820); (ii) Abell 548 is more circular near the core than 
outside the core [rj = 0.886 versus rjMKU = 0.506; and (iii) Abell 2052 has a roughly constant 
axial ratio [rj = 0.802 versus rjMKU = 0.807). 

2.5 ORIENTATION 

We also extract the cluster orientation from the Fourier expansion of the photon distri- 
bution within the annulus. 

We combine the formal and Monte-Carlo errors exactly as we do for the axial ratio. 
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Many clusters have orientations which vary as a function of radius. We combine the 
cluster orientation as a function of radius to produce an emission weighted average orientation 



where N is the number of photons within each annulus. As with and ?y, we use cluster 
simulations to test our accuracy; the third column of Figure 2 shows the test results for Og. 
The top row in the column shows the absolute error of each measurement as a function of 
the measured cluster ellipticity. As expected, the orientation error is larger for more circular 
clusters; there is no evidence of a measurement bias ((A^) = —0.85°). The bottom row 
contains a distribution of scaled errors. The best fit Gaussian is centered near zero with an 
excess width of 49% ((751 = 1. 49). Thus, we conclude that the orientation measurement is unbi- 
ased, but the uncertainties are underestimated. We increase all our orientation uncertainties 
by 49% to correct for this underestimate. 

Columns labeled Og and ao of Table 2 list the orientation angles and uncertainties for 
59 clusters. Columns labeled R['] contain the radius (in arcmin) within which we calculate 
the emission weighted average. Appendix A contains a detailed discussion of the individual 
measurements. Two clusters have variations in orientation or isophote twisting above the 3(7 



level. These two clusters are A426 {RMSg = 33.1°, (7rms = 2.4°) and A119 {RMSg = 37.0°, 



CRMS = 9.7°). In the case of an equilibrium system this twisting would be evidence of 
triaxiality; however, in these two cases there is evidence of departures from equilibrium. The 
centroid of the cooling fiow region of A426 and the centroid of the outer region of the cluster 
are significantly different (see Table 1 and MFG93). In A119, the surface brightness profile 
is inconsistent with that produced by a single mass clump. Fabricant et al. (1993) show 
that a multi-clump model suggested by the galaxy distribution produces an X-ray surface 
brightness profile consistent with the IPC image. 

Rhee & Latour (1991) analyze the X-ray orientations of 26 Abell clusters. They focus 
on the entire X-ray image, using the largest circular region centered on the peak of the 
X-ray emission which is not affected by the ribs. There are 11 measurements which appear 
in Rhee & Latour and Table 2. Although the measurement region generally differs slightly 
(Rhee & Latour typically include regions which do not satisfy our signal-to-noise constraint), 
the values generally agree; the scaled difference distribution v = {6q — Orl)/^^ has a mean 
(p) = —0.04 and a width RMSj/ = 1.2. MKU measure orientations for 49 Abell clusters. 
Because of the different region used (see previous discussion), the orientations measured by 
MKU should be consistent with our measurements only in clusters where there is no isophote 
twisting. 



There have been many attempts to understand the relationship between the radial dis- 
tribution of the X-ray emitting gas and the radial distribution of galaxies (Abramopoulos 



9o within the <(-^)> > 5 region of the cluster image. 

r n -1 





2.6 RADIAL FALL-OFF: a AND /3 
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& Ku 1983; Jones & Forman 1984). Because both the gas and the galaxies respond to the 
same gravitational potential, the Euler equation for the gas and the Jeans equation for the 
galaxies imply a relationship between the radial distributions of the two cluster components 
in an equilibrium system. If the King approximation to the isothermal sphere describes the 
galaxy distribution {pgai oc outside the core), fi describes the radial distribution of the 
gas {pgas oc r~'^^ , where (3 = ^^^2 1 T is the gas temperature, a is the galaxy velocity 
dispersion, and pnip is the mass of the average gas particle). Where the X-ray emissivity 
e oc p^^g (e-g- thermal bremsstrahlung), this /3 prescription implies a cluster X-ray surface 
brightness profile of the form in Equation 2.3. 

Some clusters in our ensemble are not well fit by this /3 model. However, the radial 
behavior of the cluster is an important morphological parameter. Thus, we examine two 
consistent approaches to measuring the radial fall-off. We measure /3 by first azimuthally 
averaging around the peak in the X-ray surface brightness to a radius of 16' (this includes 
the region within the IPC ribs for clusters located in the center of the field); we exclude 
point sources, the IPC ribs and outer regions, and separated subclumps (see Appendix A). 
We use the smoothed images to define this peak, but then produce the radial profile from 
the unsmoothed image (PSE has EWIIM~ 1.5'). The uncertainty in the surface brightness 
at a particular radius follows from the uncertainties in the image pixels which determine the 
value. We then find the parameter values which minimize the in the fit of the measured 
radial profile to a spherical /3 model. 

We require a confidence limit of 5% for the fit. In some clusters (often the ones 
with the highest quality images) this fit criterion is not met over the entire radial profile. 
This departure is sometimes caused by a cooling fiow in the cluster core; previous /3-fitting 
studies typically exclude these fiow regions (Jones & Eorman 1984). Because the X-ray image 
provides no clear boundary between the bright inner core and the outer regions, we define an 
objective core exclusion procedure. We exclude inner points in the radial profile until the fit 
to the remaining points has at least a 5% probability of consistency. EoUowing this procedure 
we include fits to all clusters except those where (i) the core radius is comparable to the IPC 
field of view, (ii) the cluster is too complex to warrant radial averaging, or (iii) we exclude 
so many inner points to get an acceptable fit that the constraints on the fit parameters are 
unacceptably weak. 

Figure 4: The accuracy of two different radial fall- 
off parameters, a and /3: From left to right the 
plots are (i) the fractional error in /3 measurements 
as a function of l3reai and (ii) the fractional error 
in the mean fall-off a versus Ureal ■ The width in 
the l3 error distribution is RM = 17% and the 
bias is {A/S/ jSreai) = —1.6%. The width in the 
a distribution is RMSa = 7.5% and the bias is 
{Aa/areal) = 0.2%. 

Here again we use the cluster simulations to test the accuracy of the /3 measurements. 
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We compare the /3 values for more than 200 pure images with the /3 values for the artificial 
IPC images (Figure 4 column 1). Clearly, the uncertainty in the measured /3's is quite high; 
although there is no compelling evidence (except perhaps for /3 > 1.2) for a measurement 
bias (( A/3//3real) = —1-6%), the scatter is quite large (RMS^ = 17%). This scatter is related 
to the correlation between the core radius Rc and the radial fall-off /3 (large Rc favors steeper 
/3— )• A3266 has a large /3 probably because of the large core radius from an ongoing merger) 
and to the limited extent of the IPC field. The Poisson noise and instrumental imperfections 
affect /3 by (i) slightly shifting the peak in the X-ray surface brightness which results in an 
increase in Rc and (ii) introducing uncertainties in the radial profile which are largest for the 
lowest surface brightness (outermost) regions. We underestimate the /3's in the clusters with 
large jS^eal in Figure 4 because the pure images are not well fit by the /3-model. Therefore, 
the cores are excluded in achieving a decent fit to the idealized radial profiles. In the artificial 
IPC images of these clusters, the Poisson noise introduces large enough uncertainties in the 
radial profiles that the entire profile (including the core) is a good fit to the /3-model in each 
case. So, in the pure images we have large Rc and /3, while in the artificial IPC images we 
have smaller Rc and /3. There is only one obvious fix for this problem: a genuinely accurate 
/3 fit requires reasonable signal-to-noise in the surface brightness profile well outside the true 
core. 

Fven with the problems associated with /3, we include measured /3's because of the 
historically attributed physical significance. (We do not include the core radii because there 
are additional problems associated with measuring them.) Columns labeled /3 of Table 2 
contain a list of 54 radial fall-off values (accurate to 17%). 48 of these values are for 
clusters with measured Lx and Tx- The KS test comparison of this sample with the Fdge 
et al. fiux limited sample shows that the two samples are statistically indistinguishable {Tx'- 
D=0.18 (40%), Lx'- D=0.20 (25%)). Figure 3 contains a histogram of the /3 values for this 
representative sample of clusters. The mean and width of this distribution are = 0.65 
and RMS^ = 0.16. The specifics of individual measurements are in Appendix A. 

Jones & Forman (1984) determine the most likely range of fi for 46 clusters observed 
with Einstein and are currently completing work on a sample of more than 300 IPC cluster 
observations (Jones & Forman 1994). 24 fi measurements appear in both the Jones & Forman 
(1984) work and Table 2, and the agreement is generally good. Using the mean of the two 
Jones & Forman /3's as the best fit value and the difference as ~ 4(7, we calculate the offset 
between our values and their values (/3 — /3jp) = —0.03. The width of the scaled difference 
distribution v = [fi — /3jp)/cr is RMSj/ = 0.75. Some of the discrepancies between the two 
measurements is undoubtedly caused by differences in our cooling-fiow exclusion strategy. 
Jones & Forman exclude inner bins until "a minimum in is reached" , while we exclude 
bins until there is at least a 5% probability that the fit and data are consistent. 

The large scatter in the /3 measurements and the fi-Rc correlation lead us to search for 
an alternative method of comparing the mean fall-off within the same physical region of 
two different clusters. We fit the mean radial fall-off a within some specified cluster region. 
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Specifically, we fit a line in log space: log(/) = a — [2a — 1) log(i?) where a is the intercept 
and a is the slope. Note that although a power law is not generally a good fit to the cluster 
radial profile (near the core), fitting a power law over a particular region (which can include 
the core) is a reasonable method of determining the mean fall-off in the surface brightness 
profile. Figure 4 contains a test of this approach using more than 200 artificial cluster images. 
The fractional error is plotted against the true value. We find no bias ((Aa/ctreal) = 0.2%) 
and the scatter in a (RMSq, = 7.5%) is much smaller than the scatter in /3. Because of the 
improved scatter, we use a in §3 to directly compare the radial behavior of cluster ensembles 
over the same physical scale. 

3. THE MORPHOLOGY-COSMOLOGY CONNECTION 

We use our ensembles of observed and simulated clusters to constrain three preferred 
cosmological models. Eor each cosmological model we ask, "Is it possible for the clusters we 
observe today to develop in a universe described by this cosmological model?" Because it 
is impossible to remove the observational effects from the observed clusters, we introduce 
these effects into the simulated data. We construct an ensemble of 24 cluster simulations (8 
sets of initial conditions evolved in three different models- the same simulations discussed 
in Evrard et al. 1993) from which we extract artificially observed ensembles similar to the 
observed cluster ensemble. 

Our observed sample consists of 52 clusters with measured X-ray temperatures and lumi- 
nosities. We exclude 5 clusters because of their low redshifts (A426, A3526, A1060, AWM7, 
and A1367). We first describe the simulations in detail and then turn to the procedure for 
comparing the simulated clusters with the observed clusters. A discussion of the results and 
the cosmological constraints follows. 

3.1 THE CLUSTER SIMULATIONS 

We employ a set of 24 simulations comprised of three sets of eight random density 

fields. The three sets correspond to three different assumptions for the underlying cosmology. 

The models we investigate are: (i) a biased, Einstein-deSitter universe with normalization 

1/2 

(Tg = 0.59 (where erg = <(((5/j//j)^)> on an Mpc scale); (ii) an unbiased (erg = 1.0), 

open universe with Oo = 0.2 and (iii) an unbiased, low density model with Oo = 0.2 and 
\o = 0.8. 

A cold dark matter power spectrum with F = ilh = 0.5 is used to generate initial 
displacement fields on a 32^ spatial grid as in Efstathiou et al. (1985). The periodic grid 
models a comoving cube of space with length L. To sample a range of cluster richness, we 
evolve two random realizations from four different box lengths L = 30, 40, 50 and 60 Mpc. 
The initial density field is constrained using the method of Bertschinger (1987) in order to 
form a significant cluster within the volume. The applied constraint is that the density field 6 
at the center of the volume, when filtered with a Gaussian spatial filter W{r) oc exp — r^/2rj 
with rf = 0.2iv, take on a present, unbiased, linearly evolved value of Sg = 3.0. This 
corresponds to perturbations 2.5 — 5 times the RMS amplitude over the range of box sizes. 
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The P3MSPH algorithm (Evrard 1988) is used to propagate the initial, linear density 
fields to their present, deeply non-linear states. Two sets of 32^ particles, coupled by gravity, 
represent the dark matter and baryonic fiuids. The baryons are approximated as an ideal, 
7 = 5/3 gas with constant mean molecular weight fi = 0.6 appropriate for a fully ionized, 
primordial plasma. The thermal energy changes adiabatically and as a result of shocks 
generated in regions of converging fiow. A cosmic baryon fraction 0;, = 0.1 is assumed for 
all the cosmologies. This approach allows for some amount of coUisionless dark matter in 
each model. 

Gravity is softened with a pairwise, Plummer potential $(r) = — l/(r^ +e^)"^/^ (in units 
where G = m = 1) with e = 0.0023iv kept fixed in the comoving frame. The smoothed 
particle hydrodynamic (SPH) method entails the use of a smoothing kernel W{r^ hi) which 
we take to be a Gaussian W{r,hi) = (vr/ij)"^/^ exp — r^//i^. The hydrodynamic smoothing 
length hi for particle i is varied to keep 



with ci = 87r. Because 95% of the kernel's weight falls within 2/ij, there are roughly 200 
particles within the domain of hydrodynamic infiuence of each particle. It is straightforward 
to show that the Gaussian kernel is similar to the B-spline employed by Hernquist & Katz 
(1989) and others with a modest rescaling of h (Evrard, Feldman & Watkins 1994). The 
smoothing length cannot increase beyond a value hmax = 0.025iv, a value derived as a fixed 
fraction of the P3M short range force cell size. An absolute lower bound on h^ of 0.5e is 
imposed. In addition, each particle's value of h^ is subject to a lower bound derived from the 
Courant condition (see Evrard 1988). In practice, the timestep is set small enough and the 
mass resolution is modest enough that this bound does not come into play. The minimum h^ 
achieved is comparable to e. The 'effective' resolution of the runs is approximately 0.005iv. 

The initial displacement fields are normalized so that a growth factor of 16 is required 
to recover the original, unbiased spectrum. Eor the biased, 0=1 models, the initial state 
corresponds to a redshift = 9. Eor the low density models, the initial state corresponds to 
redshifts = 47.5 (Oo=0.2) and 23.0 (Oo=0.2 & Ao=0.8). The code uses a fixed timestep to 
advance the particles, with 750, 1400 and 900 steps used for the 0=1, Oo=0.2 and Oo=0.2 
& Ao=0.8 runs, respectively. Energy in the Layzer-Irvine equation (Efstathiou et al. 1985) 
is conserved to typically better than 1% for the 0=1 runs and ~ 5% for the low density 
models. Because the difference equations for the gas are not accurate to second order, errors 
in the baryonic energy are larger than those of the dark matter. 



Eor each Einstein cluster observation we extract a similar observation of a simulated 
cluster from each of the cosmological models (see Eigure 5). Specifically, for each of the 
52 Einstein clusters we (i) choose one of the eight sets of initial conditions, (ii) scale the 
three (one for each cosmological model) present epoch configurations for the particular initial 




= ci 




3.2 COMPARING OBSERVED AND SIMULATED CLUSTERS 
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Emstem IPC fi=l fio=0.2 fio=0.2 & \o=0.8 





































1 

■ - 























Figure 5: Contour maps of Emstem IPC images of 4 Abell clusters and the artificial IPC images of cluster 
simulations. From the left to right the columns are (i) Emstem IPC images, (ii) artificial IPC images of 
Sl=l cluster simulations, (iii) artificial images of Slo=0.2 simulations, and (iv) artificial images of Slo=0.2 & 
Ao=0.8 simulations. From top to bottom the clusters are (i) Abell 496, (ii) Abell 644, (iii) Abell 1644, and 
(iv) Abell 2256. In producing the artificial IPC images of cluster simulations, we image the simulations at 
the same redshift as the associated Abell cluster. Contours are spaced by factors of 1.5 in surface brightness. 
For each image, the surface brightness of the bottom contour is the mean surface brightness in the outermost 
annulus which satisfies the (;^) > 5 constraint. 

conditions so that the X-ray emitting gas in the three simulations has the same temperature 
as the gas in the observed cluster, (iii) choose a random observing direction and then (iv) 
image each present epoch configuration (Evrard 1990) at the same redshift as the observed 
cluster (for Hq = 50 km s~"^Mpc~"^). The scaling conserves the mean density within the 
simulation volume while changing the temperature: T — )• subject to constant M/R^ 
implies the scale change R — )• ^/^R. For the entire ensemble of simulated images, the spatial 
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scaling varies from 0.6 to 2.9 with median values 1.1 for the 0=1 models, 1.2 for the 
Oo=0.2 models, and 1.5 for the Oo=0.2 & Ao=0.8 models. Because of the differences in 
cluster X-ray luminosity (for the default scaling, the low Oo models produce more luminous 
clusters), we select the observation time for each model so that the number of photons in the 
region of interest is identical for all four ensembles. In this way we produce four ensembles 
(1 observed and 3 artificial) of cluster observations which are homogeneous in three respects; 
each ensemble has identical distributions in (i) gas temperature, (ii) spatial scale, and (iii) 
number of cluster photons within the region of interest. This approach isolates the differences 
in cluster morphology. 
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Figure 6: Comparing X-ray Morphologies. 
From left to right each column contains 
normalized distributions of (i) the emission 
weighted centroid variation in h'^^kpc (ii) 
the emission weighted axial ratio rj, and (iii) 
the radial fall-off a. From top to bottom 
the rows correspond to the following cosmo- 
logical models: (i) fio=0.2, (iii) fio=0.2 & 
Ao=0.8, and (iii) Sl=l. The bottom row is 
the Emsiem IPC sample. Table 3 lists the 
mean and width of these distributions. 



100 300 
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To increase the sample size and fully probe the morphological characteristics of the 
simulations, we cycle through the 52 Einstein IPC observations 5 times, producing a total of 
260 different comparison-sets of cluster observations. Each set consists of one Einstein IPC 
observation and artificial observations of the same initial density field evolved in each of the 
three cosmological models (note that because we only use 52 Einstein images, each of the 
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IPC images appears 5 times). We then reduce and analyze the artificial images as though 
they were normal IPC observations. We compare the X-ray morphologies of the 4 ensembles 
of 260 images using three different parameters: the emission weighted centroid variation 
twf, the emission weighted axial ratio ?y, and the radial fall-off a. Figure 6 contains the 
distributions in these three parameters. In making these comparisons we treat each set of 
4 images as a unit to preserve the ensemble similarities in Tx, and photon number; if 
making a measurement on any one of the 4 images is impeded then we exclude the entire 
set (For example, from Table I it is clear that A2I42 has no w^^ measurement because of a 
central point source. Rather than just excluding the cluster w:^ for A2I42, we exclude all 5 
A2I42 comparison sets-4 images each.). 

3.3 THF DISTRIBUTION OF CFNTROID VARIATION 

The first column of Figure 6 contains the distribution in w:^ for the four samples of cluster 
images. For the observed clusters and the 0=1 simulations we calculate w^^ starting from an 
inner radius of 4', but for the Oo=0.2 and Oo=0.2 & Ao=0.8 simulations we exclude those 
innermost annuli with radii smaller than the scale of one "effective" resolution length. This 
is necessary because the cluster cores in the low Oq models have evolved to central densities 
so high that corresponding cooling times are less than a Hubble time. Because the present 
simulations do not include radiative cooling, the innermost cores of the low Oq simulations 
are non-physical and must be excluded. (Calculating w:^ with these bright cores included 
tends to increase the differences between the 's measured in the 0=1 model and the low 
Oo models.) Clearly, clusters with large centroid variations are much more probable in the 
observed sample and the 0=1 model than in the two low Oq models. KS test results indicate 
that the 0=1 model and the observed sample are marginally consistent, and that the two 
low Oo distributions are marginally consistent. All other combinations differ significantly. 

Because a sizable centroid variation is a signature of dynamical youth (MFG93, Mohr & 
Fvrard 1994), it is no surprise that these distributions refiect the differences in the merger 
histories of the three models. As predicted by Richstone, Loeb, & Turner (1992), the clusters 
evolving in the 0=1 model are more likely to experience significant growth at late times than 
those evolving in the low Oq models. The similarity between the 0=1 distribution and the 
observed cluster distribution suggests similarities in the evolutionary histories of these two 
ensembles: both the observed and the simulated clusters have a high probability of late-time 
growth. 

3.4 THF DISTRIBUTION OF AXIAL RATIO 

The second column of Figure 6 shows the striking differences in the distributions of i]. 
As with the centroid variations, we exclude one "effective" resolution length from the cores 
of the low Oo clusters. The absence of late-time mergers in the two low Oo models produces 
predominantly spherical clusters at the present epoch (Oo=0.2: = 0.95; Oo=0.2 & Ao=0.8: 

= 0.91). The higher probability of late-time mergers in the 0=1 models results in more 
highly fiattened clusters thus reducing the average axial ratio ((?y) = 0.70). The sample 
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of Einstein clusters = 0.80) lies intermediate between the 0=1 model and the two low 
Oo models. A KS test indicates that all four distributions differ significantly. 

The physics not included in the numerical simulations (e.g. gas cooling) will undoubtedly 
affect the distributions of rj. For example, no cooling fiows can form in the simulations, but 
observed cooling fiows are characterized by circularly symmetric peaks in the X-ray surface 
brightness. Such a fiow tends to weight the distribution of axial ratio toward 1 (especially 
for an emission weighted axial ratio). We examine this effect using X-ray images of observed 
clusters (see Figure 3). Of the 51 clusters included in the Figure 3 rj distribution, 44 have 
been analyzed for the presence of cooling fiows (Fdge, Stewart, & Fabian 1992, Stewart et 
al. 1984). 13 of these 44 clusters have central cooling times which are significantly less than 
10"^*^ years (these clusters are marked with a in Table 1 and Table 2). The dashed lines in 
Figure 3 show the rj distribution after the removal of these cooling fiow clusters. With one 
exception, the clusters containing definite cooling fiows fall in the region ?y> 0.8 (the region 
with fewer clusters in the 0=1 model). 

On the other hand, recent cluster mergers tend to weight the distribution toward the 
highly fiattened end. The simulations do follow cluster mergers properly, so including cooling 
in the simulations would shift the rj distribution toward 1. Such a shift would make the 0=1 
model more consistent with the observed clusters. Because the Oo=0.2 and Oo=0.2 & Ao=0.8 
models have rj distributions which are already more spherical than observed clusters, adding 
cooling would not significantly improve their match to observed clusters. 



Table 3. Mean (and RMS) oiw^, rj, and a Distributions 



Parameter 


Emsiem 


n=i 


fio=0.2 & Ao=0.8 


fio=0.2 


w^[kpc] 


50.1 (49.2) 


30.4 (39.3) 


6.6 (8.8) 


5.4 (7.9) 


ri 


0.80 (0.12) 


0.70 (0.17) 


0.91 (0.07) 


0.95 (0.02) 


a 


1.75 (0.32) 


1.82 (0.36) 


2.68 (0.27) 


2.88 (0.36) 



3.5 THF DISTRIBUTION OF RADIAL FALL-OFF 

The third column in Figure 6 contains the normalized a distributions for the four en- 
sembles. The fit region for a is determined separately for each of the 260 comparisons. In 
each case the inner boundary is the larger of the "effective" resolution lengths for the Oo=0.2 
and Oo=0.2 & Ao=0.8 models; the same core region is excluded from all four clusters within 
each of the 260 comparisons. As with the ?y's and 's, the radial fall-off a neatly separates 
the low Oo simulations from the 0=1 models and the observed clusters. The radial profiles 
of the low Oo clusters are significantly steeper. A KS test indicates that the a distributions 
of the 0=1 clusters and observed clusters are marginally consistent. 

It is not possible to make the low Oq clusters consistent with observed clusters by chang- 
ing the slope of the power spectrum (Crone, Fvrard, & Richstone 1994); however, gas physics 
excluded from the simulations may significantly affect the cluster radial fall-off. Metzler & 
Fvrard (1994) investigate the effects of gas ejection from galaxies in a series of simulations. 
They produce an ensemble of initial density perturbations consistent with standard, 0=1 
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Figure 7: We demonstrate the effects of ejection from galaxies on the cluster 
radial fall-off using cluster simulations which include gas ejection (Metzler & 
Evrard, 1994). The histograms contain the results of /3-fits to idealized images 
of 47 different clusters. The top histogram is the radial fall-off for those clusters 
evolved with ejection from galaxies = 0.65), and the bottom histogram 

contains the radial fall-off for the identical clusters evolved without ejection 
((/?>= 0.87). 




.4 -6 -8 1 1-2 1-4 



CDM constrained to form clusters in a range of box sizes. They then evolve each set of these 
initial conditions twice. In the first set of simulations they simply evolve the baryonic and 
dark matter components as described in §3.1. In the second set of simulations they investi- 
gate ejection from galaxies by (i) replacing 2.5(7 density peaks on galaxy scales (5 X IQ-^-^M©) 
with single, massive, dark matter particles and (ii) modeling continuous gas ejection from 
galaxies with the discreet ejection of hot gas particles from galaxies over time. They explore 
the maximal effect of ejection from galaxies by following an extreme ejection model: constant 
ejection totalling half the galaxy mass from z=^ to the present epoch. In Figure 7 we plot 
the histograms of radial fall-off fi for the clusters evolved without ejection (bottom) and 
those evolved with ejection (top). The clusters evolved with no ejection ((/3) = 0.87) have 
gas distributions which fall off more steeply than the same clusters evolved with ejection 
((/3) = 0.65); the fi distribution of the clusters evolved with ejection is very similar to the 
intrinsic fi distribution of observed clusters (see Figure 3). 

The Metzler & Fvrard simulations demonstrate that gas ejection can significantly de- 
crease the radial fall-off of the ICM density. Since the models employ a fairly extreme ejection 
history for galaxies, a more realistic case would introduce a smaller effect. Figure 6 shows 
the difference in the average value of the mean slope parameter a between the observations 
and models is ~5% for = 1 and ~40% for low Oq- Fjection at a level higher than that em- 
ployed by Metzler & Fvrard may be advocated to make the profile slopes for the low density 
runs consistent, while a more modest level of energy input from galaxies would bring the 
= 1 runs into agreement with the observations. At any rate, because the bulk of ejection 
in the low Oo models would occur at very high redshift, it would not be able to alter the 
isophote shifts or axial ratios of the present day X-ray images. The low Oq models would 
still disagree with the observations on these measures. 

3.6 THF FMISSION WFIGHTING AND CORF FXCLUSION 
To ensure that the morphological differences between the low Oq clusters and the 0=1 
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clusters are properly reflected in Figure 6, we test whether our methods of emission weighting 
and core exclusion (in the low Oq clusters) introduce measurement bias. To test emission 
weighting, we calculate and rj for the entire ensemble of Oo=0.2 & Ao=0.8 and 0=1 
simulations both with and without emission weighting. Specifically, we set = 1 in Equa- 
tion 2.2 and Equation 2.6 and recalculate the values. The comparison of more than 400 
measurements for the 0=1 and Oo=0.2 & Ao=0.8 simulations reveals that the mean change 
in axial ratio rj is 0.2% for the Oo=0.2 & Ao=0.8 clusters and 0.5% for the 0=1 clusters. 
The mean change in is 0.016 arcmin for the Oo=0.2 & Ao=0.8 clusters and 0.020 arcmin 
for the 0=1 clusters (compared to the median uncertainty for the Einstein cluster sample of 
0.083 arcmin). Clearly these changes are insignificant when compared to the large differences 
refiected in the Eigure 6 histograms. 

Measurements on the low Oq clusters are bounded on the interior by the unresolved core 
and on the exterior by the ("^) > 5 constraint. To determine whether and r] are sensitive 
to the exact size of the excluded core, we double the core and compare the measurements. 
The artificial IPC image produced from each simulated cluster is exposed to obtain the 
same number of cluster photons; therefore, doubling the diameter of the exclusion region 
increases the Oo=0.2 & Ao=0.8 observing times. Because of the increased observing times, 
the <(-^)> > 5 constraint is satisfied to larger radii. We find that for the Oo=0.2 & Ao=0.8 
simulations, increasing the size of the core exclusion region by a factor of two results in a mean 
change in of 0.026 arcmin and a mean change in rj of 0.5%. Clearly, the morphological 
measurements for this ensemble of low Oq clusters are not significantly affected by our core 
exclusion procedure. 

4. CONCLUSION 

We assemble a representative sample of 65 Einstein IPC cluster images to quantify the 
range of cluster X-ray morphologies. Eor each cluster we measure the emission weighted 
centroid variation (tWf), the emission weighted axial ratio (?y), the emission weighted ori- 
entation (^o), and the radial fall-off parameters, a and /3. We employ a Monte-Carlo 
approach to account for the effects of Poisson noise, instrumental imperfections, and fore- 
ground/background X-ray point sources. With the cluster simulations described in §3 we 
evaluate the measurement accuracy (Eigure 2 & Eigure 4). We list the measurements on each 
cluster (Table 1 & Table 2) and construct histograms (Eigure 3) of the intrinsic distributions 
in , ?y, and /3. This set of morphological measures is sensitive to the cosmological density 
parameter 0. The measurements indicate that from 50% to 70% of clusters exhibit core 
substructure, depending on the adopted confidence limit. 

We use the X-ray morphologies of observed and simulated clusters to place cosmological 
constraints. With the combined N-body/SPH code P3MSPII (Evrard 1988), we evolve eight 
different initial density fields (sampled from an ilh = 0.5 CDM spectrum) according to 3 
different underlying cosmological models: (i) 0=1, (78=0.59; (ii) Oo=0.2, (78 = 1; and (iii) 
Oo=0.2 & Ao=0.8, (78 = 1. These 24 simulations are the same ones analyzed by Evrard et al. 
(1994). EoUowing a comparison method designed to isolate the cosmological parameters, we 
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use 52 of the Einstein cluster images to compare the range of observed and simulated cluster 
morphologies. We demonstrate that (i) w^^ ?y, and a neatly separate clusters evolved in low 
Oo models from those evolved in 0=1 models (see also Evrard et al. 1994), and (ii) clusters 
like those we observe today are very rare in low Oq universes. 

How general are these results? Can any low Oq model reproduce the observed complexity 
of the X-ray images of real clusters? For the case of bottom-up Gaussian initial density fields, 
we suspect the answer is no. The reason is that merging is the most natural mechanism 
capable of generating asymmetric isophote structure in the ICM, and the present merging 
frequency of clusters is much more sensitive to the density parameter than to the shape 
of the power spectrum on cluster scales (e.g. Lacey & Cole 1993). Increasing Oq is the 
most natural way to increase the merger frequency and thereby produce a higher fraction of 
distorted X-ray images. This point is sharpened by the fact that we do not have complete 
freedom in specifying the shape of the power spectrum on cluster scales. The abundance of 
clusters as a function of X-ray temperature implies an effective spectral index n between 
— 1 and —2 (Henry & Arnaud 1991), and this range is consistent with power spectra derived 
from large-scale galaxy catalogues (Vogeley et al. 1992; Fisher et al. 1993; Feldman, Kaiser 
& Peacock 1994). 

How might the low Oq models be saved? We have seen that energy ejection from galaxies 
can fiatten the gas density profile. However, it is unclear whether a realistic level of ejection 
can produce an effect large enough to bring the a and fi measurements into agreement 
with observations. One might attempt to fix the axial ratios for the low density models 
by appealing to a higher dark matter fraction (presently 50% in our models) which could 
make the underlying cluster potentials more aspherical. However, the clusters would still be 
dynamically old, so the isophotes, even if less spherical, would retain their symmetric nature. 
Thus the centroid variation wy^ would continue to present a strong constraint. To address 
these questions directly, we intend to pursue experiments of specific low density models in 
the near future. 

Other recent attempts to constrain Oq from cluster properties also disfavor low density 
models. Matching the mean velocity dispersion of rich clusters with Oq = 0.2 requires erg ~ 
1.25 — 1.58, implying galaxies are anti-biased with respect to the total mass content of the 
universe (Lilje 1992; White, Ffstathiou & Frenk 1993). The strong increase in the fraction of 
blue galaxies in moderate redshift clusters (the "Butcher-Oemler effect") is also difficult to 
explain in a low density universe (Kauffmann 1994). A recent paper by de Theije, Katgert & 
van Kampen (1994) presents a possible exception to this trend. They find the ellipticities of 
the galaxy distributions within Abell clusters to be somewhat rounder than those found in 
their = 1 simulations. However, they did not yet have a complete set of low Oq simulations 
to compare against the observations. In general, any evidence for recent evolution or merging 
in the cluster population tends to disfavor low Oq. With improved optical and X-ray data 
combined with detailed numerical simulations, specific cases of recently merged clusters are 
becoming more apparent. Fxamples include A2256 (Fabricant, Kent, & Kurtz 1989; Briel & 
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Henry 1994), A754 (Fabricant et al. 1986, Zabludoff & Zaritsky 1994) and Coma (Burns et 
al. 1994). 
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APPENDIX A 

In this appendix we list the details associated with measurements on specific Emstem IPC cluster im- 
ages. Reductions for those clusters not listed follow the standard procedure described in the text. We define 
a set of guidelines and attempt to apply them uniformly to all cluster images. In measuring the cluster 
radial fall-off, these guidelines include (i) excluding well separated subclumps from the radial profile where 
possible, (ii) avoiding radial averaging for those clusters which are very complex, (iii) excluding point sources 
from the radial profile, and (iv) truncating the radial profile well before the IPC ribs. For these guidelines 
are (i) excluding regions with bright point sources (the fainter, less obvious point sources are accounted for 
with a Monte-Carlo procedure), (ii) including subclumps except in those cases where there is independent 
evidence that the subclump is a chance superposition and not dynamically related to the main cluster, and 
(iii) avoiding those clusters where the ^-^^ > 5 region is so small that fewer than 4 centroids are measured. 
For rj and Og we (i) avoid regions containing bright point sources and (ii) exclude well separated subclumps. 
Angular distances and directions are given with respect to the peak in the X-ray surface brightness. 

A85: a//?: subclump to S excluded. r\jdo'- Truncated at i? = 5.4' to exclude subclump. subclump included. An 
extensive radial velocity study (Malumuth et al. 1992) provides evidence that 7 galaxies are associated with a 
foreground group; however, the spatial distribution of these 7 galaxies is not well centered on the subclump in 
the X-ray emission. Therefore, we consider the subclump to be evidence of an ongoing merger and therefore 
include it in the measurement. 

A133: a Iii: three point sources excluded (18' N, 13' SW, and 18' S). »?/6io/«)3: (f ) > 5 region does not include 
these three point sources. 

A168: ajji: cluster too complex for radial analysis. rilOo/wg: no rib or point source contamination within (■^) > 5 
region. 

A400: a//?: point source 15' S of cluster excluded. No /? value included because core radius from /9-fit larger than 

IPC field of view. rj/Oolw^: point source 15' S of cluster is outside (■§■) > 5 region. 
A426: 'qjBol'UJx- truncated at i? = 13.9' to avoid rib contamination. 

A478: a//?: excluded two point sources (lO' NW and 25' NE). rj/Oolw^: no point sources within (■§■) > 5 region. 
A539: a//?: excluded one point source (3 W/NW). 'qjBol'UJx- point source inhibits measurements. 
A592: ajfi: truncated at i? = 13.3'. wj: (-f) > 5 region too small for measurement. 

A754: ajji: point source excluded (25 SW). No ji included because core radius is comparable to IPC field of view. 
rilOo/wg: no point source within (■^) > 5 region. 

A780: a If): point source excluded (17' S/SE). rj/Oo: point source outside (■^) > 5 region, wg: (■§■) > 5 region too 

small to make measurement. 
A1060: rilOo/w^: truncated at i? = 13.9' to avoid rib contamination. 

A1367: a//?: IPC field of view inadequate to measure /?. r^ldolw^: point source near peak in surface brightness 

removed. Truncated at i? = 12.8' to avoid rib contamination. 
A1644: a//?: excluded subclump to NE. rijBo: truncated at i? = 6.4' to exclude subclump. w^: subclump included. 

An extensive radial velocity study failed to detect a superposed subclump; therefore we consider the subclump 

to be evidence of an ongoing merger. 
A1650: ajji: one point source excluded (20' SE). wr- {§) > 5 region too small for measurement. 
A1656: rj/Oo/wg: Analysis truncated at i? = 14.9 to avoid rib contamination. 
A1689: wr- {§) > 5 region too small for measurement. 

A1736: ct/lS: excluded two point sources (17' SE and 30' E). Tj/do/wg: point sources are exterior to the (■^) > 5 
region. 

A1775: ct/lS: excluded three point sources (12' SE, 22' SE, and 25' S/SE). J]/6o: point source external to (■^) > 5 
region, wg: (■§■) > 5 region too small for measurement. 
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A1795: a//?: excluded point source (5' SW). rijdojws'- point source inhibits accurate measurements. 
A1983: wr- {§) > 5 region too small for accurate measurement. 

A2029: ct/f): excluded two point sources (3' SE and 15' S/SW). Tj/do/wg: point source inhibits accurate measure- 
ments. 

A2052: a//?: excluded one point source (15' NW). rj/Oo/wg: point source is exterior to (■^) > 5 region. 
A2063: ct/f): excluded one point source (3' S/SE). Tj/do/wg: point source inhibits accurate measurement. 
A2065: a//?: excluded two point sources (10' NW, 15' W). rj/Oo/wg: point sources are outside (■^) > 5 region. 
A2124: a//?: excluded possible point source (13' NW). rj/Oo/wg: point source outside (■^) > 5 region, wg: (■^) > 5 

region too small for measurement. 
A2142: ct/lS: excluded two point sources (4' NE, 17' SE). Tj/do/wg: point source inhibits accurate measurement. 
A2147: a//?: excluded three possible point sources (20' E/SE, 15' S/SE, 20' NW). No /? included because no fit with 

at least a 5% probability of being consistent with the data exists. rj/Oo/wg: point sources lie outside (■^) > 5 

region. 

A2151: a//?: system too complex for radial averaging. rj/Oo- truncated at i? = 6.4' to exclude subclumps E of main 

clump, wg: included subclumps. 
A2199: ct/lS: excluded one point source (18' NW). rj/Oo/wg: point source external to (■^) > 5 region. 
A2255: a//?: excluded two point sources (13' N, 20' NW). iildolwg: point source are outside > 5 region. 

A2256: /?: no /9-fit with at least a 5% probability of being consistent with the data exists. 

A2410: a//?: too complex for radial averaging, r^: complex nature made measurement impossible, wg: both clumps 

but no point sources included in ^-^^ > 5 region. 
A2440: a//?: system too complex for radial averaging, r^: complex nature made measurement impossible, wg: 

^■^^ > 5 region too small for measurement. 
A2593: a/f): excluded two point sources (7' N, 15' SW). rj/Oo'- point sources outside (■^) > 5 region, wg: (■§■) > 5 

region too small for measurement. 
A2597: wg: {§) > 5 region too small for measurement. 

A2626: a//?: excluded one point source (7' NE). wg: (■^) > 5 region too small for measurement. 

A2634: a//?: excluded 2 point sources (15' N/NW, 23' S/SW) and subclump (12' NW). Radial velocity and pho- 
tometry study provides evidence that subclump is due to background cluster (Pinkney et al. 1993). r^ldo/wg: 
truncated at i? = 8.5' to exclude background cluster and point source. 

A2657: a//?: excluded two point sources (7' W, 15' SW). ri/Oo/wg: truncated at i? = 6.4' to exclude point sources. 

A2670: a//?: excluded four point sources (10' N, 15' N, 15' E/SE, 20' E). »?/6io: point sources are outside (f ) > 5 

region, wg: (■§■) > 5 region too small for measurement. 
A2877: a//?: excluded point source (2'SW). ri/Oo/wg: bright point source near center of cluster makes accurate 
measurements impossible. 

A3158: a/f): excluded two possible point sources (15' S, 12' E/SE). rj/Oo/wg: point sources are outside (■^) > 5 
region. 

A3186: Wg: (■§■) > 5 too small for measurement. 
A3266: a//?: large core radius contributes to large /?. 

A3376: ct/f): excluded one point sources (20' S/SW). rijBojwg: no rib or point source contamination within (■^) > 5 
region. 

A3391: a//?: excluded one point source (13' N/NE). rj/Oo/wg: point source outside (■^) > 5 region. 

A3395: a//?: bimodal structure makes radial averaging pointless. rijBo'- considered each subclump (NE and SW) 

separately, wg: Peak of the brighter (SW) subclump used as the true cluster center. 
A3526: rijBojwg: truncated at i? = 9.6 to avoid rib contamination. 

A3532: ajji: excluded extended emission near ribs to W. wg: (-^^ > 5 region too small for measurement. 
A3667: ajfi: excluded three point sources (12' NW, 20' E/SE, 15^ S/SW). ^jdoj-wg: truncated at i? = 7.5' to avoid 
point source contamination. 

CygA: a//?: excluded subclump 12' NW. r\jdo'- truncated at i? = 8.5' to avoid subclump. wg: subclump included. 
MKW3S: wg: {§) > 5 region too small for measurement. 

3C129: ce/P: excluded one point source (15' N). No /? because the best fit core radius is comparable to the IPC field 
of view. J]/6o: truncated at i? = 11.7' to avoid point source contamination. 
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